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Constructing appropriate unitary matrix operators for new quantum algorithms and finding the 
minimum cost gate sequences for the implementation of these unitary operators is of fundamental 
importance in the field of quantum information and quantum computation. Evolution of quantum 
circuits faces two major challenges: complex and huge search space and the high costs of simulating 
quantum circuits on classical computers. Here, we use the group leaders optimization algorithm 
to decompose a given unitary matrix into a proper-minimum cost quantum gate sequence. Using 
this procedure, we present the circuit designs for the simulation of the Toffoli gate, the amplifica- 
tion step of the Grover search algorithm, the quantum Fourier transform, the sender part of the 
quantum teleportation and the unitary propagator of the Hamiltonian for the hydrogen molecule. 
The approach is general and can be applied to generate the sequence of quantum gates for larger 
molecular systems. 



I. INTRODUCTION 

Quantum computation promises to solve fun- 
damental, yet otherwise intractable problems in 
many different fields. It is commonly believed 
that advancement in quantum computing sci- 
ence will bring new polynomial time algorithms, 
and some NP-complete problems shall be solv- 
able in polynomial time, too. To advance the 
quantum computing field, finding circuit de- 
signs which can execute algorithms on quan- 
tum computers (in the circuit design model of 
quantum computing) is important. Therefore, 
it is of fundamental importance to develop new 
methods with which to overcome the difficulty 
in forming a unitary matrix describing the al- 
gorithm (or the part of the computation), and 
the difficulty to decompose this matrix into the 
known quantum gates Realizing the theo- 
retical problems of quantum computers requires 
the overcoming decoherence problem [2\. Re- 
cently, West et al. demonstrate numerically 
that high fidelity quantum gates are possible 
in the frame work of quantum dynamic and de- 
coupling 0]. 

The problem in the decomposition of a given 
unitary matrix into a sequence of quantum logic 
gates can be presented as an optimization prob- 
lem. Williams and Gray [4| suggested the 



use of genetic programming technique to find 
new circuit designs for known algorithms, and 
also presented results for quantum teleporta- 
tion. Yabuki, Iba 5] and Peng et al. Q fo- 
cused on circuit designs for the quantum tele- 
portation by using different genetic algorithm 
techniques. Spector [7[ explains the use of ge- 
netic programming to explore new quantum al- 
gorithms. Stadelhofer [8( used genetic algo- 
rithms to evolve black box quantum algorithms. 
There are also some other works [11-QjJ which 
evolve quantum algorithms or circuits by using 
genetic programming or genetic algorithms. Re- 
view of these procedures can be found in [l2j . 

Evolution of quantum circuits faces two ma- 
jor challenges: complex and huge search space 
and the high costs of simulating quantum cir- 
cuits on classical computers. Previous works 
focused on specific test cases and small scale 
problems. In this paper, we use our recently 
developed group leaders optimization algorithm 
(GLOA)[13[-an evolutionary algorithm-to de- 
compose the unitary matrices representing some 
quantum algorithms or the unitary propagator 
of a given many-body Hamiltonian. Using two 
construction methods to find the unitary matrix 
representation of control gates, we have been 
able to efficiently use various kinds of quantum 
gates in the optimization process without in- 
creasing the computation time. In addition to 



explaining the details of the approach, we show 
the result circuit designs for the operators of the 
Grover search algorithm; the sender part of the 
quantum teleportation; the Toffoli gate; and the 
quantum Fourier transform (while some of the 
results are new circuit designs, some of them are 
known circuit designs). For the simulation of 
the hydrogen Hamiltonian, we give two circuit 
designs which have much less cost than the cir- 
cuit design found by mapping the Hamiltonian 
from fermionic operators to the Pauli operators 
& 

This paper is organized as follows: after defin- 
ing the objective function to be minimized, we 
describe the algorithm used in the optimization; 
sec II VI are devoted the implementation details; 
and in Sec [V] the result circuit designs are pre- 
sented and discussed. 



II. DEFINING OBJECTIVE FUNCTION 

The definition of the objective function is an 
important factor in the optimization process. In 
our case, there are two factors which need to 
be optimized within the objective function: the 
correctness and the cost of the circuit. Let Uf 
be the matrix found by the optimization and U g 
is the given unitary matrix to be decomposed. 
The correctness is defined in fijl as: 
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where C is the correctness; N = 2" [n is the 
number of qubits); the symbol f represents the 
complex conjugate transpose of a matrix; and 
Tr(..) is the trace of a matrix. This definition 
ignores the global phase differences which are 
physically indistinguishable; hence, makes the 
optimization easier by diversifying the reach- 
able solutions for a problem in complex space. 
C obtained from Eq.([T]) determines how much 
Uf looks like U g and is in the range [0, 1] since 
the product of two unitary matrices is another 
unitary matrix all eigenvalues of which have ab- 
solute value 1, and when Uf = U g , \Tr(UfU^)\ 



is equal to 2"; hence, the value of the correct- 
ness is 1. 

The cost of a circuit describes the level of ease 
with which this circuit is implemented; in order 
to make the implementation of a circuit easier 
and the circuits less error-prone, the cost of a 
circuit also needs to be optimized by minimizing 
the number of gates in the circuits. Hence, in 
addition to the correctness, the objective func- 
tion should include the cost of the circuit. 

However, defining the cost of a circuit is not 
an easy task due to the fact that each quan- 
tum computer model may have a different cost 
for a given quantum gate. Here, we simply 
take into consideration the number of qubits 
on which a gate is operating and the number 
of control qubits the gate has; these are the 
common factors increasing the implementation 
cost of a gate in all quantum computer mod- 
els. Therefore, because quantum operations on 
distant qubits may not be feasible to imple- 
ment physically [161 ] , for the implementations of 
the circuits which requires the control of many 
qubits it may be possible to separate the circuit 
(the whole computation) into smaller parts for 
the partial implementation; in this sense, hav- 
ing quantum gates operating on closer or adja- 
cent qubits in the design eases the implementa- 
tion difficulty. 

Considering the reasons in the previous para- 
graph, we define the cost of different types of 
quantum gates as follows: The cost of a single 
gate is 1 and less than the cost of any type of 
the control gates since the implementation of 
a single gate on quantum computers is much 
simpler than the implementation of the control 
gates. To estimate the cost of the control gate, 
instead of placing a constraint on the target and 
control qubits by requiring them to be adjacent 
as done in we take the cost of a gate with 
the closer target and control qubits; the number 
of qubits between the control and target qubits 
is multiplied by a constant. This constant is 2 
for the regular control gates and 3 for the multi- 
control gates. For instance, while the cost of the 
CNOT gate whose control and target qubits are 
neighbors is 2 x 1 = 2 (2 is the constant and 1 
is the number of qubits between the target and 
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control qubits.), the cost for the Toffoli gate is 
3 x 2 = 6 (3 is the constant for the multi-control 
gates, and 2 is the number of qubits between the 
first control and the target qubits.). After find- 
ing the cost of each gate in a circuit, the cost of 
the whole circuit is defined as the sum up of the 
costs of the gates in the circuit. The following 
expression finds the cost of the circuit in FigfTa} 

Cost = 2x1 + 2x1 + 2x2 + 2x1 = 12. (2) 

The general objective function is defined by 
including both the correctness and the cost of 
the circuit with some weights: 
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where the constants a and ft are the weights for 
the correctness and the cost, and defined as: 



0< (a,ft) < 1; 
a + ft = l. 



(4) 



The cost in the objective function is also scaled 
to the range [0, 1] by taking ft less than 1. The 
expression aC must be dominant in the objec- 
tive function in order to find exact or more ac- 
curate solutions. Because of the definition of 
the cost and the correctness, a should be much 
greater than ft to make aC dominant. There- 
fore, the choices of a and ft have effects on the 
number of iterations which is necessary to find a 
desired solution for a problem. Taking a = 0.9 
and ft = 0.1 is the optimum choice which re- 
duces the number of iteration and increases the 
correctness by making the correctness dominant 
in the objective function and giving higher pri- 
ority to the correctness than the cost in the op- 
timization process. Hence, in the experiments, 
we set a = 0.9 and ft — 0.1. 

It is important to note that the value of the 
objective function never becomes zero since all 
circuit designs with at least one gate have a cost 
value. For the test cases, in addition to the val- 
ues of objective functions, we give the C values 
for the circuits to make the cost and correctness 
values in the objective function more coherent. 
Also note that on the results of the objective 
function there may be some computer round-off 
errors. 



III. GROUP LEADERS OPTIMIZATION 
ALGORITHM 

Group leaders optimization algorithm 
(GLOA) described in more detail in [l3[ is 
a simple and effective global optimization 
algorithm that models the influence of leaders 
in social groups as an optimization tool. The 
general structure of the algorithm is made 
up by dividing the population into several 
disjunct groups each of which has its leader 
(the best member of the group) and members. 
The algorithm which is different from the 
earlier evolutionary algorithms and the pivot 
method algortihm jl7H19 | consists of two 
parts. In the first part, the member itself-the 
group leader with possible random part-and a 
new-created random solution are used to form 
a new member. If the formed new member is 
a better solution to the problem than the old 
member, it replaces the old one. This mutation 
is defined as: 

new member — r\ "portion of old member 
U r2 portion of leader 
U r3 portion of random, 

(5) 

where n , r%, and r% determines the rates of the 
portions of the old member, the group leader, 
and the new-created random solution; which 
form the new member, and sum to 1. 

In addition to the mutation, in each iter- 
ation for each group of the population one- 
way-crossover (also called the parameter trans- 
fer) is done between a chosen random member 
from the group and a random member from 
a different-random group. This operation is 
mainly replacing some random part of a mem- 
ber with the equivalent part of a random mem- 
ber from a different group. If new formed mem- 
bers give better solution to the problem, then 
they survive and replace the old members of the 
groups; otherwise, they do not. The amount of 
the transfer operation for each group is defined 
by a parameter called transfer rate. The val- 
ues of the parameters of the algorithm used in 
the experiments are given in the next chapter. 
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The flow chart of the algorithm for our opti- 
mization problem is drawn in FiglHJ The more 
information about the algorithm can be found 
in ref.[I3. 



IV. IMPLEMENTATION DETAILS 
A. Representation of Quantum Gates 

In the optimization, we use similar represen- 
tation method to the method of Cartesian Ge- 
netic Programming |2(| in which each function 
set and inputs (in our case, gates and qubits) 
are represented as integers and genotypes in- 
cluding the inputs and the gates are represented 
as integer strings. The difference is: since the 
quantum gates may have an effect on the whole 
system, the gates should be represented in time 
steps; that means we cannot give the same in- 
puts to two different gates at the same time as 
it is done in classical circuits. In the string of 
the genotypes each four numbers represent the 
gate, the qubit on which gate operates, the con- 
trol qubit, and the angle for the rotation gates. 
The integers for the gates are determined by 
looking at the index of the gates in the gate set. 
For a gate set {V, V*,Y}, an example numeric 
string representing the circuit in Fig llal is as: 

2 3 2 0; 3 2 1 0; 2 3 1 0; 1 3 2 0; 3 2 1 0, 
where the each group of four numbers separated 
by semicolons describes the each quantum gate 
in the circuit: the first numbers with bold fonts 
identify the gates, the last number is the an- 
gle, and the middle integers are the target and 
the control qubits, respectively, (the semicolons 
and the bold fonts do not appear in the real 
implementation) . 



B. Construction of the unitary matrices 
for the gates 

Each quantum gate is considered as a unitary 
matrix which multiplies an input unitary ma- 
trix (in our cases an identity matrix): UiUi npu t- 
Hence the sequence of Ui matrices represent the 



sequence of gates. For each evaluation of the ob- 
jective function, a unitary matrix of order N for 
each gate in the circuit need to be constructed 
by taking the Kronecker product of the unitary 
matrix Ucate and the required identity matri- 
ces. To find the matrix representation of dif- 
ferent type of control gates, Ucate, we use two 
similar pseudo codes which are run in O(N) and 
O(l) on a given identity matrix of order N and 
help to reduce the computation time of the op- 
timization. 

The main idea of these pseudo codes is: to 
use the order of qubits in the state register rep- 
resented as \qoqi.-.q n ) (n is the number of qubits 
on which the control gate acts); find the related 
indices for these qubits in a given identity ma- 
trix; and using the distance between the con- 
trol and target qubits in the state register, to 
change the related elements of this matrix with 
the elements of the controlled elementary gate. 
The initial positions of the control and target 
qubits in the identity matrix are located by us- 
ing the order of the control and target qubits 
in the state register, \qoqi---q n )- And then, the 
jumping amount from this initial positions in 
the unitary matrix are determined by using the 
distance between the control and target qubits, 
and the elements in the related positions are 
changed according to the elements of the single 
elementary gate that acts on the target qubit, 
and is controlled by the control qubit. 



1. Pseudo code to construct the unitary matrix of 
a single- control quantum gate 



Input(target, control, U) 

{Either the control or the target qubit is zero.} 
{U is a 2 d+1 by 2 d+1 identity matrix.} 

^. -i- ^target. 

2: c = t + 2 contro1 ; 
3: m = 2 4 - 1 - 1; 
4: for k = to m do 
5: i = c + 2k; 
6: j = t+ 2k; 

7: Uu = u 00 ; Uij = uoi; U Jt = u w ; Ujj = u u ; 
{where Moo, woii an d Uu are the 
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matrix elements of the single gate.} 
8: end for 

return U ; 

The above pseudo code builds a unitary matrix 
which represents the different types of the con- 
trol gates possessing only one control and one 
target qubits. It locates the related elements of 
the identity matrix, and changes them with the 
elements of the elementary gate acting upon the 
target qubit. 

In the code, d- which is the absolute difference 
between the control and target qubits-helps to 
find the number of qubits covered by the control 
gate and is the key element in the determina- 
tion of the number of times the for-loop shall 
run. The matrix U is initially defined as a 2 d+1 
by 2 d+1 identity matrix which shall be the rep- 
resentation of the control gate at the end of the 
algorithm. The variables c and t are the indices 
used to find the i and j which determine the in- 
dices of elements to be changed inside the loop. 
After finding the correct indices, the four ele- 
ments of U are changed in each iteration with 
the elements (moo, woij w io, an d wn) of the sin- 
gle gate that operates on the target qubit. Since 
the loop runs m + 1 times, the total 4 x (m + 1) 
elements are replaced by the elements of the sin- 
gle gate. At the end of the algorithm, the ma- 
trix U represents the desired unitary matrix for 
the control gate, and is given as an output of 
the algorithm. 

The running time of this pseudo code is O(N) 
in the worst case since the loop, which is the 
dominant in the running time, runs rn+l times. 



2. Pseudo code to construct the unitary matrix of 
a multiple- control gate 



Input {target, control, U) 

{Either the control or the target qubit is zero.} 
{U is a 2 d+1 by 2 d+1 identity matrix.} 

1: d = \control — target\; 

2- £ ^target. 

3: c=t + 2 contr ° l ; 



4: TO = 2 d ~ 1 - 1; 
5: i = c + 2m; 
6: j = t + 2m; 

7: Uu = uqq; Uij = uqi', Uji = itio; Ujj = u\\\ 
{where Moo, uoi, u io, and Un are the 
matrix elements of the single gate.} 

return U; 

Unlike in the regular-control gates, all of the 
qubits between the control and target qubits 
are also the control qubits in the multi-control 
gates. The idea of this pseudo code is almost 
the same as the previous pseudo code, but the 
absence of a for-loop makes this simpler than 
the first one and run in 0(1). 



V. RESULTS AND DISCUSSIONS 

In the experiments, some parameter adjust- 
ments have been done as follows: The default 
gate set consists of the rotation gates, the Pauli 
operators, the square root of not (V), the com- 
plex conjugate of V gate (V*), and the con- 
trolled type of these single gates. For the matrix 
representation of single quantum gates used in 
this paper, please refer to 0, HH, [22[. While in 
the case of the unitary propagator of the Hamil- 
tonian of the hydrogen the angle is defined as a 
multiple of 0.005 and in the range [0, 2ir], in the 
other cases it is defined as kir, and k is multiple 
of 0.125 and in the range [0, 2]. During the eval- 
uation of the objective function in the optimiza- 
tion, all quantum gates, except single gates, are 
formed by using the pseudo codes given in the 
previous section. 

The parameters for the optimization algo- 
rithm are given in Table H The correctness (C) 
and the minimized objective function values and 
the number of iterations are noted in Table [TT] 
with respect to all circuit design results that are 
given in the paper. 

The circuit diagrams in figures are drawn by 
using the output of the optimization program. 
As an example, the output of the program for 
the decomposition of the unitary matrix for the 
diffusion part of the Grover search algorithm is 
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as follows: 
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where the column G represents the name of the 
gate, T is the target qubit, C is the control 
qubit, and Q determines the angle values for 
rotation gates. The number of rows is equal to 
the maximum number of gates. The zeros are 
the control gates with the equal control and tar- 
get qubits, or a field the gate does not need to 
use such as the control qubit field for the single 
gates. Fig|2al shows the circuit diagram for this 
output. 

The test cases are divided into two categories: 
the known quantum algorithms (the details of 
these algorithms can be found in and the 
unitary propagator of the hydrogen Hamilto- 
nian. In the former case, the result circuit de- 
sign does not show much significant improve- 
ment over the known circuits; therefore some of 
them have the same length or the structure as 
the known circuit designs since a great deal of 
work has been done on these algorithms, and 
these cases except three qubit quantum Fourier 
transform can be considered as simple cases for 
the optimization because of the small size of the 
circuits which are formed by only non-rotation 
gates (the angles for rotation gates are also pa- 
rameters to be found by the optimization algo- 
rithm; which makes the optimization procedure 
harder.). However, implementation difficulties 
of different circuit designs with the same cost 
and length may not be the same in different 
quantum computer models; hence, finding dif- 
ferent circuits is important for the implemen- 
tation. In addition, the results show the effi- 
ciency and the ability of the optimization pro- 
cedure (see Table |n| . The result circuit designs 
for the quantum algorithms are shown in FigQ] 
FigH FigEJ and FiggJ The case of the hydro- 
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FIG. 1: The found circuit designs (a), (b), and 
(c) for the Toffoli gate. 



gen Hamiltonian which finds the low cost circuit 
designs is explained in detail in the following 
subsection. 



The Hamiltonian of the hydrogen 
molecule 



It has been shown that the ground and ex- 
cited state energies of small molecules can be 
carried out on a quantum computer simula- 
tor using a recursive phase-estimation algorithm 
[23M25| . Lanyon et al. reported the application 
of photonic quantum computer technology to 
calculate properties of the hydrogen molecule 
in a minimal basis [l4j . Here, we show how our 
method can be used to reduce the number of 
gates needed to perform the simulations. 

Fermion model of quantum computation is 
defined through the spinless fermionic annihi- 
lation (flj) and creation (aj) operators for each 
qubit j (j=l, ... , n), where the algebra of 2n 
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FIG. 2: The found circuit designs (a) and (b) for the two-qubit amplification part of Grover 

search algorithm. 
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FIG. 3: The found circuit designs (a) and (b) for the two-qubit quantum Fourier transform. 



elements obey the fermionic anti-commutation 
rules [H: 



{ai,aj} = 0, {a,i,a]} = 6ij, 



(7) 



where {A, B} = AB + BA defines the anti- 
commutator. Using the Jordan- Wigner trans- 
formation, the fermion operators is mapped to 
the standard quantum computation operators 
through the Pauli spin operators [26|: 



\fc=l 



(8) 



Once the electronic Hamiltonian is defined in 
second quantized form, the state space can eas- 
ily be mapped to qubits. The electronic Hamil- 
tonian in second quantized form is described as 



m tea tea: 



H = h pq a p a q + - ^ h pqrs a p a q a r a s , (9) 



pqrs 



where the integrals h pq and h pqrs evaluated dur- 
ing the Hartree-Fock procedure are defined as 
follows: 



h pq — 



Xq(x) 

(10) 



and 



dx\dx2 



X p (x 1 )x* q (x 2 )Xr(x2)Xs(xi) 



ri2 



where 



(11) 

is the distance between the a th nu- 
cleus and the electron, r\2 is the distance be- 
tween electrons, V 2 is the Laplacian of the elec- 
tron spatial coordinates, and Xp( x ) is a selected 
single-particle basis. Whitfield et al.[27| consid- 
ered H in Eq.© as + . Since h pqrs = 

foqprsi and SO hijji — hjiij — ^ijij — ^jijii 

they noted the parts of the Hamiltonian as fol- 
lows: 



ffW = hna\ai + h22a\ a '2 + ^3303^3 + haa\a4, 

(12) 

and 
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FIG. 4: The found circuit design for the three-qubit quantum Fourier transform 
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FIG. 5: The found circuit designs for the sender part of the quantum telportation. 
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(13) 



Using their findings [13] for the spatial inte- 
gral values for atomic distance 1.401a. it., the 
Hamiltonian and the unitary propagator for this 
Hamiltonian, e~ lHt (t was taken as 1), are found 
as matrices of order 16. The algorithm is run 
with the parameter values in Table U and Table 
HIl for the unitary propagator, and the decom- 
posed circuit design is shown in Figl6] It is im- 
portant to note that the circuit in Figj6]does not 
include the approximation error coming from 
the Trotter-Suzuki decomposition; however, it 
has some small errors which can be gaged from 
the values of the correctness and the objective 
function in Table HIl 

In addition to the exact unitary propaga- 
tor, the unitary operator at bond distance 
1.401 atomic units is found by simulating the 
approximated-general circuit designs for the 
simulation of the Hamiltonian in the paper 
[13] (In 0, the non-commuting terms in the 



Hamiltonian are approximated by following the 
Trotter-Suzuki decomposition, and the given 
approximated circuit design is a general circuit 
for the unitary propagator; the values of the 
angles in the quantum gates depend on the spa- 
tial integral values.). The decomposition of this 
unitary propagator is found as in Fig with 
the value of correctness 0.9961. The circuit de- 
sign in FigEl consists of 7 quantum gates while 
the complete circuit design in [27j consists of 87 
quantum gates. 



VI. CONCLUSION 

To be able to simulate Hamiltonians of atomic 
and molecular systems and also apply quantum 
algorithms to solve different kinds of problems 
on quantum computers, it is necessary to find 
implementable quantum circuit designs includ- 
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ing the minimum cost and number of quantum 
gate sequences. Since deterministic-efficient 
quantum circuit design methodology is an open 
problem, we applied stochastic evolutionary op- 
timization algorithm, GLOA, to search quan- 
tum circuit designs for given unitary matrices 
representing algorithms or the unitary propa- 
gator of a molecular Hamiltonian. In this pa- 
per, in addition to explaining the ways of the 
implementation and design of the optimization 
problem, we give circuit designs for the Grover 
search algorithm, the Toffoli gate, the quantum 
Fourier transform, and the quantum telepor- 
tation. Moreover, we give the two circuit de- 
signs for the simulation of the Hamiltonian of 
the hydrogen molecule by decomposing the uni- 
tary matrix operators found by following the 
fermionic model of quantum computation and 
simulating the circuits given in [27|. In the 
case of the hydrogen molecule we reduced signif- 
icantly the number of gates needed to simulate 
the unitary operator from 87 quantum gates to 
7 quantum gates. The approach is general and 
can be applied to generate the sequence of quan- 
tum gates for larger molecular systems. 
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FIG. 6: The circuit design for the unitary propagator of the Hamiltonian of hydrogen molecule. 
The unitary propagator is found by using the spatial integral values in [27| and the definitions for 
the annihilation and creation operators in Eq.® into the Eq.©, Eq. (fT2j) . and Eq. (fT3|) . 
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FIG. 7: The circuit design for the unitary propagator of the hydrogen Hamiltonian. The unitary 
propagator is found by simulating the circuits and pseudo-code given in [27] . 

TABLE I: The values of the parameters within the group leaders optimization algorithm 



Parameter 
name 


Number of groups 


Number of popula- 
tion in each group 


ri 


T2 


^3 


Number of transfer 
for each group 


Value 


15 


25 


0.8 


0.1 


0.1 


Number of Variables* i 


2 



fSince each gate has four variables, the total number of variables is equal to four times maximum number of gates. 



TABLE II: Optimization parameters and results for the test problems 



The unitary matrix 


Number of 


Maximum 


Number of 


Objective 


Correctness 




qubits 


number of 


iteration 


function 


(C) 






gates 




value 




Toffoli Gate, Figllal 


3 


8 


500 


0.09167 


1 


Toffoli Gate, Figflbl 


3 


8 


500 


0.09286 


1 


Toffoli Gate, Fig[lc] 


3 


8 


500 


0.09286 


1 


Sender part of teleportation, Figl5al 


3 


8 


2500 


0.08752 


1 


Sender part of teleportation, Figl5bl 


3 


8 


2500 


0.08752 


1 


Diffusion Operator of GSA, Fig|2al 


2 


8 


500 


0.08335 


1 


Diffusion Operator of GS, Figllbl 


2 


8 


500 


0.08571 


1 


Two-qubit QFT, Fig|3al 


2 


8 


500 


0.08752 


1 


Two-qubit QFT, Figlgbl 


2 


8 


500 


0.08752 


1 


Three-qubit QFT, Figg] 


3 


12 


2500 


0.09565 


1 


Simulation of the H2 Hamiltonian, Fig[6] 


4 


12 


10000 


0.1081 


0.9870 


Simulation of the H2 Hamiltonian, Fig[7] 


4 


12 


10000 


0.0892 


0.9961 
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Create n number of groups, generate and 
evaluate initial population for each group. 






1 > 


f 




Choose the leaders for each group from the 
group members. 


> 


r 



Update all members by using Eq, (18). 

(A random solution also needs to he created for Eq. (18), for the 
update of each member) 



A given-number of times, do parameter transfer 
(One-way crossover) for each group from other 
groups . 




FIG. 8: The flow chart of the group leaders optimization algorithm 
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